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ABSTRACT 

Recent observations have provided strong evidence for the presence of an Electron 
Scattering Region (ESR) within the central regions of AGNs. This is responsible for 
reprocessing emission from the accretion disk into polarised radiation. The geometry 
of this scattering region is, however, poorly constrained. In this paper, we consider 
the influence of gravitational microlensing on polarised emission from the ESR in 
the quadruply imaged quasar, Q2237+0305, demonstrating how correlated features 
in the resultant light curve variations can determine both the size and orientation 
of the scattering region. This signal is due to differential magnification between per- 
pendicularly polarised views of the ESR, and is clearest for a small ESR width and 
a large ESR radius. Cross- and auto-correlation measures appear to be independent 
of lens image shear and convergence parameters, making it ideal to investigate ESR 
features. As with many microlensing experiments, the time-scale for variability, being 
of order decades to centuries, is impractically long. However, with a polarization fil- 
ter oriented appropriately with respect to the path that the quasar takes across the 
caustic structure, the ESR diameter and radius can be estimated from the auto- and 
cross-correlation of polarized light curves on much shorter time-scales. 
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1 INTRODUCTION 

While quasars are amongst the most luminous sources in 
the Universe, being at cosmological distances ensures that 
their small angular size renders them unresolved to mod- 
ern telescopes. While direct imaging of the central engines 
of quasars is not possible, more novel techniques are able 
to reveal the scales of structure of several prominent fea- 
tures [e.g. reverberation mapping c an reveal the si ze of the 
enveloping Broad-Line Region; see IPetersonl (| 1998b ]. Grav- 
itational microlensing is one such method, occurring when 
the light from a distant source is magnified by the pres- 
ence of stellar masses crossing the line of sight. Within the 
Galaxy, where single or binary stars pass in front of more 
distant sources, microlensing has been used to great effect 
in studying the surf ace properties of dista nt stars, such as 
limb darkening (e.g. iRhie &; Bennettlll999h . 

Microlensing also occurs in a cosmological context, 
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where the light from a distant quasar has been multi- 
ply imaged by a foreground galaxy. As with the Galac- 
tic case, stars within the lensing galaxy which cross the 
line of sight to these images can induce substantial mag- 
nification effects. Unlike the Galactic case, the number of 
stars influencing a particular image can be many thou- 
sands, and hence the resultan t pattern of magnification 
can be extremely complex (see Kayser, Ref sdah fe StabelJ 
1986; IWambsganss, Paczvnski. &: Katdll99(j for examples). 
Regions within the quasar will respond differently to the 
magnifying stars, dependent upon their size and loca- 
tion, and hence the light curve variations as the stars 
cross the line of sight encode the underlying quasar struc- 
ture. Understanding and, more importantly, deconvolving 
this encoding has been the focus of a substantial amount 
of study in recent years, and has included the develop- 
ment of computational approaches to account for the grav- 
itational lensing by a large population of sources (see 
ISchneider. Kochanek, & Wambsganssll2006l for a recent re- 
view). Microlensing has been used to study structure on 
a range of scales within quas ars, both limiting the size 
of the central accretion disk (|Yoneharal l200ll : | Bate et al.l 
l2008l : lBlackburne et alj|2010h and including the properties of 
the "big blue bump" ([Wambsganss, Paczvnski, Sz Schneiderl 
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Figure 1. Schematic model of an AGN viewed from the side, 
with the jet-emitting core and accretion disk surrounded by an 
electron scattering region. Beyond this is a broad line emission 
region that may overlap the ESR and is itself surrounded by a 
dusty torus. In the top left corner is a top- view schematic that 
depicts how unpolarised light from the central source becomes 
polarised after scattering by the ESR; the arrows represent the 
direction of E polarisation. 



Il990h . On larger scales, microlensing of the b road line 
regio n has been the subject of t h eoretical (iNemirofj 
1988; Schnei der fc Wambsganssl Il990l : lAbaias et aTn"2 002: 



Lewis & Ibat al l20Q4f) and observa t ional ^Richards et al.l 
20041 : IWavth. O'Dowd. Websterl l2005h investigations, 



with differential microlensing between the inner disk and ex- 
tensive line em itting region provid ing further clues to quasar 
structure fe.g. iKeeton et al1l2QQ6h . 

In this paper, we will consider the influence of mi- 
crolensing on the prominent electron scattering region 
(ESR) responsible for the repro cessing of quasar accret ion 
disk emission into polarised light ([Kishimoto et al.l 2008a, b). 
Although the use of polarimetry in probing the Broad Ab- 
sorption region through microlensing ha s been suggested 
(jBelle fc Lewisll2000l : [Hales fc Lewisll2007h . here we consider 
the influence of microlensing on the scattering region in gen- 
eral. Evidence for the microlensing of an electron scattering 
region in the gra vitationally l ensed quasar H1413+1143 has 
been provided bv lChae et al.1 (|200lh , where it was found that 
one of the images had a rotated polarisation position angle 
and a somewhat higher degree of polarisation than the other 
images; this was attributed to microlensing of an ESR. 

The focus of this paper is to examine the influence of 
gravitational microlensing on fiducial models of the ESR, 
using numerically generated maps of the complex magnifi- 
cation patterns found in quasar microlensing. The structure 
of this paper is as follows; we detail our model of an electron 
scattering region in Section [2] as well as describe the math- 
ematical formalism behind gravitational microlensing and 
the method by which we apply this to our model. Section [3] 
presents the results of cumulative and time-series statistical 
analyses with the aim of discovering methods of identify- 
ing the size and orientation of an ESR. The conclusions are 
presented in Section [4] 



2 BACKGROUND AND APPROACH 
2.1 Source Model 

Quasars represent the most powerful of active galactic nu- 
clei, with emission arising from very energetic compact 
regions thought to be supermassive black holes powered 



by accretion disks. Although there are several models for 
the specific aspects of quasar structure (e.g. lElvisI [2000; 
ISchild fe Vakulikll2003h . the consensus is that a broad line 
emitting region (BLR) lies outside the accretion disk and a 
dusty torus exists beyond that, both of which fuel the cen- 
tral accretion disk. The dusty torus is typically considered 
to be on the scale of 0.1 kiloparsecs, but sizes for both the 
BLR and torus are not well constrained. Recent publications 
about the dusty torus have even argued for a min imum ra- 
dius on the order of parsecs (e.g. lAgol et~ai1l2009h . 

It has been suggested that an additional electron scat- 
tering region (E SR) is required to explain polarized emission 
from q uasars ([Taniguchi &; Anabukil [l999 ; IKishimoto et al.l 
2008a) as polarimetry appears to cut through the scattered 
emission of the dusty torus, implying that scattering from 
the interior is able to produce spectra of the accretion disk 
without contamination from dust emission. However, the ge- 
ometry of the ESR an d its position relative to t he BLR is still 
a subject of debate ([Kishimoto et al1l200 8b); a combined 
model possessing all of the key features within a quasar is 
displayed in Figure [1] 

The quasar's accretion disk is represented by a 2D 
Gaussian of radius 



r 1/2 = 4 x 10 14 (A/zim) 1 - 5 m, 



(1) 



where 7*1/2 is the half-light radius of the accretion disk and 
A is the emission wavelength in the quasar's rest frame 
(Ago l et al.ll2009h . For the purpose of this investigation, 
A = 0.7 fim (i.e. infrared) is used t o comply with the ESR 's 
scattering wavelengths of interest ([Kishimoto et al1l2008bh . 
The accretion disk's presence is primarily to place a lower 
bound on the ESR radius. 

For the purpose of this study, we ignore the BLR as, 
although po larisation levels from this region are potentially 
quite high ([Antonuccil 1 19991 ) , it should be possible to dis- 
count its ef fects due to the presence of emission line spectral 
signatures ([Kishimoto et al.ll2Q08ah . The dusty torus is also 
neglected in our model due to its immense size, meaning 
that it will be largely unaffected by micro lensing, producing 
inten sity variations of less than 1% (e.g. iRefsdal & Stabell 
Il99lh . In contrast, the ESR is considered to be contained 
within distances below the minimum radius of the dusty 
torus and is thus far more susceptible to small-scale magni- 
fication variations that affect its intensity. 

As a fiducial model, we assume the case where the 
quasar jet is aimed at the observer and the ESR is seen as 
a rotationally symmetric, radially Gaussian annulus; while 
simple, this model captures the key physical features of the 
ESR, although we intend to generalize this model in fu- 
ture contributions. With this assumption, the source is an 
extreme-case Type 1 object, where the nuclear region inte- 
rior to the dusty torus is visible. The polarisation E vector 
of scattered emission is tangent to the circumference of the 
annulus (as in the insert, Figure [1} and, in the case of a 
face-on disk, any circularly symmetric ESR will result in 
zero net polarisation. However, due to the effects of differ- 
ential magnification, parts of the polarised emission will be 
preferentially magnified, resulting in a net polarisation. A 
polarising filter may then be used to view the annulus such 
that either the horizontal or vertical component of emission 
is detected; these will be called the "horizontally polarised" 
(HP) and "vertically polarised" (VP) filters. The filter may 
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Figure 2. An example of a simple ESR model with minimum 
radius of 0.64 pc and a width of 0.64 pc. The top left image is the 
ESR viewed through a HP filter rotated anti-clockwise by 7r/6, 
i.e. Of = 7r/6 in Equation [2] If the polarisation filter was not 
rotated, i.e. if Of =0, the wings would lie horizontally. The top 
right image is the ESR viewed through a VP filter rotated anti- 
clockwise by 7r/6. Without a rotated filter, the wings would lie 
vertically. Red and blue in these two panels represent high and 
low intensities, respectively. The bottom left image is the total 
flux of the ESR. The bottom right image is the total polarisation 
of the ESR as defined in Equation [3] Red and blue in this panel 
represent high and low levels of PolTot? respectively. Axis units 
are in Einstein radii for Q2237+0305 (see Section [272]) . 



then be rotated, admitting light polarised at an angle from 
the horizontal or vertical. These points are summarised in 
the following equation, describing the visible ESR intensity 
profile I as a function of polarisation filter and orientation, 
annulus dimensions, annulus intensity and location over 
the source plane: 



I = Io sin 



if) exp 



(r-(r E + w E /2)Y 
2(^/4) 2 



(2) 



where r is the radial distance from the quasar centre, r E is 
the minimum radius of the ESR, we is the width of the ESR, 
Io is the flux normalisation, is the source plane location, 
6 a indicates use of a horizontal or vertical filter, and Of is 
the rotation of the filter. The third factor in the equation 
implements the Gaussian cross-section of the annulus, which 
peaks at the average radius of the ESR and emits 95% of its 
total photon flux from within its "width" . The second factor 
implements the effect of a polarising filter. A horizontal or 
vertical filter is specified by a = or ty/2 respectively, and 
rotation of the filter is specified with a non-zero value oiOf. 
In all cases the visible source profile is a pair of mirror-image 
"wings" , as depicted in Figure [2] The top left image is the 
ESR viewed through the HP rotated filter; the wings would 
be horizontal except for the fact that the filter is rotated by 
angle Of = 7r/6, leading to a corresponding profile rotation. 



The top right image uses a rotated VP filter; the wings would 
be vertical except for the fact that the filter is at Of = 
7r/6. In these top two panels, red and blue indicate high and 
low intensities, respectively. The orientation of the profile is 
therefore specified by one of HP or VP plus the angle Of. The 
bottom left image is the total flux profile for the annulus, 
and the bottom right image is the total polarisation, denned 
as 

/hp — I VP 



PoIto 



(3) 



/hp + /vp 

at Of = 7r/6, where /hp and /vp are the intensities of the 
horizontally and vertically polarised components, respec- 
tively. In this panel, red and blue indicate high and low 
values of Polxot , respectively. 



2.2 Gravitational Lensing 

The formali sm of microlensing is well descri bed in the lit- 
erature (see ISchneiderT Ehlers. fe Falcol lT992): however, we 
reproduce the key aspects here. Strong gravitational lensing 
occurs when the light from a distant source is deflected by 
a massive object, allowing several light paths to connect it 
to an observer and resulting in multiple images of the dis- 
tant source. For galaxy-mass lenses, the typical separation 
of these images are of order an arcsecond, although individ- 
ual stars within the lensing galaxy can also induce additional 
deflections as they pass through the line-of-sight to a source. 
While these deflections, of order a micro- arcsecond, are too 
small to produce additional resolvable images, they can pro- 
duce substantial brightness variations; this is the observable 
consequence of gravitational microlensing. 

For this study, the key characteristic length that arises 
in gravitational lensing is the Einstein Radius; this describes 
the radius of a ring seen around a point-mass lens when there 
is perfect alignment between the source, lens and observer. 
Projected onto the source plane, this is given by 



Co 



4GM D os D, s 



D , 



(4) 



where M is the mass of the lens and D xy refers to the angular 
diameter distance between x and y; the subscripts s, /, and 
o representing source, lens, and observer respectively. 

While gravitational lenses are three-dimensional mass 
distributions, the cosmological distances involved mean that 
in practice we can treat them with the thin lens approxima- 
tion, where the mass distribution is projected onto a lens 
plane, giving a surface mass density £(£) at point £ in the 
plane. This is normalised by a critical surface mass density 
of 

c 2 D os 

such that the normalised surface mass density is given by 
a = S/S cr , also referred to as convergence (a). The mapping 
of intensity between the source and image planes can be 
expressed with a Jacobian matrix of the form 



A - 



■ (J — 71 -72 
—72 1 — a + 71 



(6) 



where 7 = \J ^\ + 7! is called the shear and incorporates 
the long-range influence of the lensing mass distribution. 
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0.47 
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0.38 


0.43 
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0.65 


0.68 


(-)2.94 
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0.59 


0.56 
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Table 1. Convergence, shear and mean theoretical magnification 
for the four images seen in the gravi tational lens Q2237+0305 
(take n from the best-fit model of iRix. Schneider. &; Bahcalll 
Il992h . The sign of the magnification denotes parity. 
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Figure 3. Two different examples of magnification maps for im- 
age B with a = 0.38 and 7 = 0.43. The second one is rotated by 
<j> = — 7r/4. Spatial units are in Einstein Radii for a solar mass 
star. 



The theoretical magnification of an image is then given 



fith = (det A)~ 



(1 - a) 2 - 7 2 



(7) 



Note that formally infinite magnifications are possible when 
the denominator of Equation [7] goes to zero. This defines a 
network of critical curves in the image plane that map to 
caustics in the s ource plane ; patterns of light and dark seen 
in Figure [31 ( see IWittlll990l . for examples), described in full 
later in the text. 

For convenience it is common to set 72 to zero, thus 
choosing an orientation where the dominant caustic struc- 
tures are considered to be horizontal (as an examp le, see Fig- 
ure 2 in IWambsganss, Paczvnski, &; Katzlll990h . However, 
rotation of the caustic network is possible by noting 



tan(20) = 

7i 



(8) 



where 4> represents the angle between the dominant caus- 
tic structures and the horizontal axis (Figure [3]). We will 
henceforth refer to these structures as caustic "bands" . 

In practice, each image of a multiply-imaged lensed 
quasar is modelled with a convergence and shear value de- 
scribing the image magnification within the overall model 
of the system. As the lens and source move relative to each 
other over time, microlensing by the stellar objects within 
the lens will deviate the magnification above and below this 
value. For the purposes of this study, we will focus upon 
the gravitational lens Q2237+0305 (iHuchra et all Il985h 
due to its significant hist orical monitoring (|Wozniak etHrtl 
l200d : lUdalski et alJl20o3 ). which has revealed ongoing mi- 
crolensing. Q2237+0305 is an attractive laboratory for mi- 
crolensing investigations because of the low redshift of 
the lensing galaxy, lead ing to high microlensing variability 
(|Kavser &; Ref sdal 1989). It is a quadruply imaged z — 1.69 
quasar seen through the bulge of a z = 0.04 spiral galaxy; 
the Einstein Radius for a solar mass star in this system is 
~ 0.06 pc. The adopted values of the shear and convergence 
for each image in Q2237+0305 are given in Table [U 

Analysis of recent Q2237+0305 observations by the 
Spitzer Space Telescope used several m odels to constrai n 
the minimal radius of the dusty torus (|Agol et alJ 1 20091 ) . 
The largest value w as calculated as 3.83 pc for an interstel- 
lar medium model (|Draindl2003h . so we treat this as the 
conservative upper bound on the ESR radius. 



2.3 Numerical Analysis of Gravitational 
Microlensing 

For this study, we us ed ba ckwards ray -t racin 
dKavser. Refsdal Stabefl Il98& IWambsganssl fl99S„ 
iGarsden fc Lewis! l2010h to generate microlensing magnifica- 
tions; here, rays are fired from an observer through a plane 
of lensing masses. These masses deflect the path of the 
rays, which are then traced into a distant source plane. The 
plane is represented as a pixel grid, with the number of rays 
that land in each pixel of the source plane corresponding 
to the magnification of a pixel-size source situated at that 
pixel. The result is a magnification map (Figure [3} of 
the source plane, where bright areas indicate locations 
where the source will be highly magnified and dark areas 
indicate where they will be demagnified. Extended sources 
are represented as a flux profile pixel grid and convolved 
with the magnification map to produce a map for the 
microlensing of that source. 

The analysis in this paper considers all microlensing 
masses to be 1 Solar mass, with the Einstein Radius (ER) 
projected into the source plane of 0.06 pc. Note that, as 
the Einstein Radius is proportional to y/~M, all distances 
related to this problem are scalable by this factor. While 
any realistic field of microlensing bodies is expected to pos- 
sess a mass distribution, the overall statistical properties 
are proportional to the square root of the mean mass (e.g. 
iLewis fe Irwi3ll995h . 

All magnification maps are 150 ER (9 pc) across, with 
the inner 75 ER (4.5 pc) used so as to avoid edge effects. 
Considering that the ESR is e xpected to lie ~ 0.1 pc awa y 
from the centre of the quasar (jTaniguchi fe Anabukilfl999h . 
we henceforth set the inner radius of the ESR to be 0.08 
pc (1.33 ER) and the annulus width to be 0.32 pc (5.33 
ER), unless otherwise stated. This model will be referred to 
as the "base" or fiducial ESR. To obtain a statistical sam- 
ple, we generate 20 maps for each of the four lensed images 
in Q2237+0305 (see Table H), each map having the lensing 
masses at different randomly selected positions. Every map 
is also rotated by <j) — 0, — 7r/12, —tt/6 and — 7r/4, as enabled 
by Equation [8] 
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3 RESULTS 

In the first subsection we present the results of magnifica- 
tion map analysis for different source profiles. The second 
subsection studies the dynamical effects of microlensing by 
analysing light curves. 



3.1 Cumulative Statistics 

Magnification probability distributions of magnification 
maps are shown in FigureUfor the base ESR as seen in image 
B, using unrotated maps and as viewed through HP and VP 
filters at two filter orientations each; unrotated, and rotated 
by Of = 7r/4. The distributions appear very approximately 
Gaussian in shape, centred on the expected theoretical mag- 
nification. In Figure [4al the top and bottom panels use un- 
rotated HP and VP filters, respectively. Figure [4bl presents 
the same data but for Of — tt/4. By subtracting the former 
distribution from the latter, Figure l4cl shows that a relative 
angle of 7r/4 between polarisation filter and the band- like 
caustic structures, evident in Figure [3] biases the distribu- 
tion towards the mean with less extremes in magnification. 
However, subtracting the distribution in Figure I4bl from a 
histogram for Of — tt/2 (not shown) generates a vertical in- 
version of Figure [4c] with an increased proportion of high 
and low magnification with respect to the average. Further 
simulations confirm this trend for any additional rotation of 
tt/2. 

The explanation of this requires consideration of how 
the HP and VP wings overlap with the structure of the mag- 
nification map. If the relative orientation is such that either 
the HP or VP wings are parallel to the caustic bands, there 
is a high probability that either the wing pair crosses caus- 
tics and magnification is larger than average, or the wings 
are both off the caustic structure and magnification is below 
average. Presumably, a perpendicular orientation also leads 
to such extremes; if one wing lies across a caustic then the 
other is likely to cross the same band, and likewise if one 
wing is off the structure then the other is too. In contrast, 
a relative angle of ty/4 does not double the magnification 
effect for a wing pair via simultaneous traversal of the same 
caustic, nor can the pair completely overlap a band as they 
do when parallel to the caustic structure. Hence the prob- 
ability of magnification is more tightly centred around an 
average. 

Figure [5] extends this argument by displaying the joint 
probability distributions of the base ESR VP filter mag- 
nifications (y-axis) and HP filter magnifications (x-axis) 
for Q2237+0305, all as seen in image B at varying ro- 
tations of the polarisation filter and magnification maps. 
Joint distributions have been fruitfully used i n other works 
dLewis fc Ibatal I2QQ4I : iBrewer fc Lewisl 120051 : lAbaias et all 
l2007h and indicate the probability of HP and VP filtered 
magnifications occurring at the same source location, i.e. 
map point, and thus how well magnifications vary in "lock- 
step". They therefore indicate expectations of one magnifi- 
cation based on the other, and can constrain models from 
observations. The probability value is indicated by colour, 
with blue being and red "most likely" . All the panels rep- 
resent the same measure, but for different orientations of 
the magnification map (0) and the various rotations applied 
simultaneously to both HP and VP filters (Of). 
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Figure 4. Magnification probability distributions (histograms) 
for Q2237+0305 B using the base ESR, unrotated maps, (a) rep- 
resents the use of unrotated HP and VP filters, (b) represents a 
rotation of Of = 7r/4. (c) shows the additional number of counts 
that the histograms in (b) have over (a). 
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Figure 5. Joint distributions for the magnifications of the base ESR seen through the VP filter (vertical axes) versus the magnifications 
seen through the HP filter (horizontal axes), all as seen in image B with varying map and filter rotations. Different columns correspond 
to various magnification map rotations and different rows correspond to various polarisation filter rotations. White lines denote equal 
magnification of both components. 



All distributions appear symmetric across the lines of 
equal magnification. However, the top left and bottom right 
plots show widest distributions when the relative angle be- 
tween filter and magnification map is a multiple of 7r/2. The 
spread of the distribution is a result of the extreme relative 
magnifications that occur when one wing pair of the ESR 
lands on a caustic and the other pair does not, be it parallel 
or perpendicular. In contrast to the and ty/2 relative angle 
plots, the four on the diagonal from bottom left to top right 
represent a relative angle of ty/4 between the VP wings and 
the shear axis. They show the greatest level of correlation 
between the magnification of horizontally and vertically po- 
larised light as both wing pairs are symmetric across the 
axis perpendicular to the caustic bands; there is no position 
for the ESR with respect to the magnification map that al- 
lows significant overlap between the caustic structure and 
one wing pair but not the other. 

Figure [6] replaces the variation in cj> (here is every- 
where) with the different lensed images. It shows that there 
is a tightening of the distribution around the line of equal 
magnification for relative angles of 7r/4 and this is a general 
effect for any of the four images. As expected, the distribu- 
tions are also centred around the theoretical magnifications 
for each lens, with values presented in Table [1] 

Thus far we have only commented on features that in- 
dicate the relative orientation between the shear axis of a 
lens and a polarisation filter. However, it is also clear that 
source size plays a role in the probability distribution of 
magnification. Figure [3 shows distributions for image D with 



0f = <)> = and different ESR annulus sizes, indicated by 
r and w. These sizes indicate t hat a small ESR on the low 
end of the assumed size range (jTaniguchi fc Anabukilll999l ) 
has a long tail in its magnification distribution, suggesting 
large fluctuations in relative intensity. A heightened degree 
of correlation between observed horizontal and vertical po- 
larisation intensities should be expected due to the small 
source size and the large amount of time the entire annulus 
is on or off a caustic relative to the amount of time a po- 
larised component is. In contrast, a wide ESR annulus only 
ever has minor portions overlapping caustic bands, and thus 
variations in the intensity around the theoretical mean are 
much smaller. 

For annuli of large inner radius, Figure [3 shows that the 
distributions are tightened around a mean. The explanation 
for the change is similarly to do with the ESR area; the 
magnification map has a fixed characteristic length scale for 
caustics, as implied by Figure [3l and thus a larger ESR will 
have an intensity that is more robust to perturbations in 
magnification as it moves relative to the lens. 

We finally note that the last panel in Figure [3 is an anal- 
ogy to the dusty torus situation that we excluded from the 
model in Section \2. II There is a lack of correlation between 
the VP and HP regions due to their size, and magnification 
variability is on the order of only 5%. Hence, if the intensity 
of observed polarisation does not fluctuate significantly, it 
implies the existence of a very large ESR. 

With the spatial separation between HP and VP wings 
of the model ESR annulus, as evident in Figure [2] there is an 
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Figure 6. Joint distributions for the magnifications of the base ESR viewed through the VP filter (vertical axes) versus the magnifications 
view through the HP filter (horizontal axes), with unrotated maps and for images A, B, C, D. Different columns correspond to various 
polarisation filter rotations and each row corresponds to a different image. White lines denote equal magnification of both components. 



expectation that the caustic structure of magnification maps 
would preference either the magnification of horizontally or 
vertically polarised light. The symmetry of distributions in 
Figs. [5] and \6\ suggests that no such effect is discernible for 
a default inner radius of 0.08 pc and a width of 0.32 pc. 
However, the low ESR- width panels of Figured demonstrate 
that there is a clear asymmetric preferencing effect if the 
ESR annulus is sufficiently thin. 

Figure [8] investigates this further by fixing an ESR 
width of 0.04 pc and varying the inner radius. The left col- 
umn has Of = so the HP wings are parallel to the caustic 
bands. The right column has Of = tt/2 so the HP wings 
are perpendicular to the bands. It is evident in the bottom 
two panels that sufficiently large inner radii allow greatest 
fluctuations in magnification for the parallel components. 
This is because the perpendicular wings are less influenced 
by perturbations from traversing individual caustic bands; 
they lie across a spread of caustics and only have a rela- 
tively small portion of their area crossing new bands as the 
lens moves. The significance of this is that the shape of the 
joint distributions gathered over a substantial period of time 
will potentially describe both the width of the ESR and its 
distance from the accretion disk. 



3.2 Time-Series Statistics 

In Q2237+0305 the estimated transverse ve- 
locity of the lensing galaxy is 600 km s _1 



([Wambsganss, Paczvnski, fe Schneider] Il990h : leading 
to a projecte d velocity across the source plan e of ~ 6500 
km s _1 (see iKavser, Refsdal, &: StabelJ Il986h . It is this 
large projected velocity, and correspondingly small vari- 
ability time-scale, that makes Q2237+0305 such an ideal 
laboratory in which to study gravitational microlensing. 
Using this, time-series statistics that predict periodicities 
can be tested and the structure of a quasar's ESR may be 
resolved by comparing simulated results with observations. 

The primary tool for temporal statistics in microlensing 
is the light curve, as seen in Figure [9] The left column dis- 
plays five horizontal light curves for image A, (f> = 6f = 0, 
from a single map. The blue line represents the source profile 
for the total annulus flux, the green line relates to the HP 
filtered intensity, and the red line represents the VP filtered 
intensity. The right column is the same, except for verti- 
cal paths. By taking a cut through a generated convolution 
map, we represent the observed intensity of a source as its 
origin moves along the relevant path of the magnification 
map. It is reasonable to expect from the results gathered 
in Section I3.ll and Figure [9] that fluctuations due to caustic 
crossings may be periodic in nature if the caustic networks 
themselves are spatially correlated. 

We employ simple cross- and auto-correlation mea- 
sures to examine whether periodicities occur in the light 
curves, which would indicate an approximately regular in- 
teraction between the caustic str ucture and both the size 
and orientation of the ESR fsee ISeitz &; Schneider! fl994; 
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Figure 7. Joint distributions for the magnifications using the VP 
filter (vertical axes) versus the magnifications using the HP filter 
(horizontal axes), all as seen in image D with unrotated maps and 
unrotated filters. The size of the ESR model changes from panel 
to panel; the rows correspond to the changing inner radius of the 
ESR annulus and the columns to its width. White lines denote 
equal magnification of both components. 



ISeitz, Wambsganss, &: Schneider! [l994l for an in-depth dis- 
cussion of microlensing correlation functions). Figure [10] 
shows auto-correlations of light curves for the VP filter of 
image B only, using horizontal map paths (left panel) and 
vertical map paths (right panel), and with cj) = Of = 0. The 
source here is a small ESR with inner radius = 0.02 pc and 
width = 0.04 pc. Horizontal light curve paths are parallel 
to the caustic bands and this is reflected in the left panel 
which shows a smooth auto-correlation. Vertical paths are 
perpendicular to the caustic bands and show a less smooth 
auto-correlation. As the effect of relative rotation between 
shear axis and polarising filter merely reduces the degree 
of correlation, we will only display correlations for vertical 
light curves in what follows. Upon examination of the rele- 
vant correlation functions for the other images A, C and D 
(not shown) , there appears to be repetitive structure within 
the light curves in steps of approximately 5.83 ER (0.35 
pc). Statistical analysis over more samples will be required 
to confirm this result. 

Figure [TT] shows auto- and cross-correlations for verti- 
cal light curves from all images of the base ESR, using total 
flux and both the HP and VP filters, with = f = 0. The 
rows represent different images and the columns different 
correlation measures. The meanings of the symbols are: "F" 
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Figure 8. Joint distributions for the magnifications of the VP 
filter (vertical axes) versus the magnifications of the HP profile 
(horizontal axes), for image D with unrotated maps. The rota- 
tion of the polarising filter and the ESR radius vary, with the 
ESR width set to 0.04 pc. The columns correspond to relative 
orientations of the polarisation filter and the rows correspond to 
the inner radius of the ESR annulus. White lines denote equal 
magnification of both components. 



as total ESR flux (annulus), "HP" as the HP filter, "VP" 
as the VP filter, "Auto" as auto-correlation, and "X" as 
cross-correlation. These panels reveal that each of the four 
lens images generate similar correlations for the base ESR, 
indicating that convergence and shear values in this range 
do not affect the light curves in any significant way. Thus 
the regions of correlation are likely to be directly related to 
the shape of the model annulus. More importantly, the size 
of the source is estimable from the correlation curves, with 
the auto-correlation showing a distinct peak after a lag of 
approximately 71 years (0.48 pc). This value is also the dis- 
tance between the edge of the annulus that first encounters 
a caustic and the inner edge of the back half that follows, 
suggesting that variations in magnification due to one point 
will be observed again when the source's back edge crosses 
the point. 

The cross-correlation curves of Figure [11] also give un- 
usual peaks at approximately 327 years (2.2 pc), which may 
imply periodic structure in the caustic networks as the dis- 
tance scale is larger than that pertaining to the source. More 
obvious though is the sharp peak correlating horizontally 
and vertically polarised magnification at distances of about 
33 years (0.22 pc), corresponding roughly to the vertical 
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Figure 11. Correlation curves for vertical light curve cuts for the base ESR, with unrotated maps and polarisation filter, for all images. 
From left to right, the columns display the magnification auto-correlation of total flux (annulus), HP filter, VP filter, and the cross- 
correlation between total flux and HP filter, total flux and VP filter, and between HP and VP filters. All possible vertical light curves in 
all unrotated maps were used, and the results averaged. "F" is total ESR flux (annulus), "HP" pertains to the HP filter, "VP" pertains 
to the VP filter, "Auto" indicates auto-correlation, and "X" indicates cross-correlation. The vertical axis is in arbitrary units and the 
horizontal axis is in years; for auto-correlation this represents the lag between light curves. 



distance between the average position of the VP wings and 
either tip of the HP pair. 

To further investigate whether the correlation features 
corresponding to ESR size are general, Figure I12al shows 
correlation measures and Figure [l2bl displays a diagram in- 
dicating the ESR separations that relate to the lags for the 
cross- and auto-correlations in Figure [T2al The left panel in 
Figure [T2al shows the HP filter light curve auto-correlation 
and the right panel shows the cross-correlation of HP filter 
light curves with VP filter light curves; all with <j) = Of = 0. 
"Auto" indicates auto-correlation, and "X" indicates cross- 
correlation. Overlaying each plot is either a dashed line that 
corresponds to the time the ESR takes to travel its diam- 



eter (2te + we) or a dotted line that relates to the radius 
(rs + we/2), both defined graphically in Figure fl2bl There 
is good agreement between both the lines and the corre- 
lation peaks; if one HP wing crosses a caustic, it is highly 
probable the second will too once the ESR travels the length 
of its diameter. Likewise, the perpendicular VP wings will 
also cross the caustic once the ESR travels the length of its 
diameter. 

The curvature of the wings at larger inner radii and 
widths means that the correlation length is lower than the 
ESR diameter, shown by the weakened correlation between 
dashed lines and peaks in the bottom three panels of Figure 
I12al More detailed calculations relating to the spatial mean 
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Figure 9. Light curves for the base ESR and image A, using a sin- 
gle unrotated map and unrotated polarisation filter. The left panel 
shows five light curves from horizontal paths, the right shows five 
vertical paths. The paths were separated by 1/4 of the map width. 
The blue line represents the source profile for the total annulus 
flux, the green line for the HP filter, and the red line for the VP 
filter. Distance units are in Einstein Radii. 
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Figure 10. Auto-correlation curves using the VP filter, with un- 
rotated maps, for image B. A small ESR of 0.02 pc for inner 
radius and 0.04 pc for width is used as the source. All horizontal 
or vertical paths are used, and the results averaged. The left panel 
shows the results for horizontal light curves and the right panel 
displays the vertical analogue. The horizontal axis is in Einstein 
Radii and represents the lag used for auto-correlation measures. 
The vertical axis is in arbitrary units. 



of wing intensity will improve the agreement. However, it is 
evident that the size of an ESR according to our model is 
discernible from the relevant correlation curves. 




(a) 



) 



(b) 

Figure 12. Correlation measures using vertical light curves, with 
unrotated maps and polarisation filter, for image D. The ESR 
model varies. In (a) the left panel is an auto-correlation using the 
HP filter, the right panel is a cross-correlation between HP and 
VP filters. "Auto" indicates auto-correlation, and "X" indicates 
cross-correlation. Dashed lines mark lags corresponding to the 
ESR diameter and dotted lines denote lags corresponding to ESR 
radius (te + we/2 according to Section fe.lj) . The vertical axis 
is in arbitrary units and horizontal axis is in years, (b) displays 
the ESR distances that relate to the lags for the cross- and auto- 
correlation in (a). 



4 CONCLUSION 

In this paper we have presented a numerical study of the in- 
fluence of gravitational microlensing upon the electron scat- 
tering region thought to exist in the inner regions of quasars. 
Focusing upon the quadruply imaged quasar, Q2237+0305, 
we considered a fiducial model with a 0.08 pc inner ra- 
dius and a 0.32 pc width con sistent with recent observations 
(jTaniguchi fe Anabukill 19991 ). As seen in Figure [5] the mag- 
nification of the two polarised annulus components can also 
be well correlated if the relative angle between filter and 
caustics is an odd multiple of 7r/4. It was also discovered 
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that, by orienting a polarising filter and defining an axis 
of polarisation as parallel to the magnification map's shear 
axis, larger ESR dimensions allow for an increased level of 
events where one polarisation is magnified more than the 
other. 

By considering the cross- and auto-correlation of the 
resulting light curves, it is seen that periodicities in the mi- 
crolensing magnification maps are imprinted on the light 
curve, providing further clues to the ESR scattering ge- 
ometry, with the size of the ESR being discernible from 
the caustic-crossing auto-correlation curve of the polarised 
"wings" that are parallel to the caustics, or the cross- 
correlation between the magnification of both polarisations. 

Polarimetric observations during a known microlens- 
ing event would provide more information than constant 
polarimetric monitoring, which is observationally more ex- 
pensive. Q 2 237+030 5js regularly monitored photometri- 
cally ([Udalski et al.l 120061 . e.g. the OGLE project) and a 
photometric microlensing event could "trigger" an intense 
polarimetric study [such over rides have bee n proposed for 
Q2237+0305 in the past (e.g. I Webster] l2004h ]. which would 
directly probe the structure of the ESR. 

It must be noted that, to obtain long term variability 
statistics, photometric and subsequent polarimetric moni- 
toring will be required on decade to century time-scales, 
even to observe the characteristic correlation lags indicated 
in Figure I12al for the small-scale fiducial model. Neverthe- 
less, a monitoring campaign on that order of duration would 
already be able to compile a magnification range indicative 
of ESR size, as suggested by Figure [71 especially for a small 
inner radius. Constraining a large-scale ESR should be pos- 
sible too, as the preferencing effects for polarizing filter ro- 
tations, as shown in Figure are likely to be observed over 
the same order of time. 

For this study, we considered a fiducial physically mo- 
tivated model of the ESR, represented as a circular torus 
about the central accretion disk. While we have examined 
some of the freedom of this model, there is substantial scope 
to extend considerations to more complex geometries (such 
as a tilting of the emission regions relative to the line of 
sight). We will consider these in further contributions. 
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